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16 Abstract 
ane 

‘Engineering data compatable with ‘the {computer program ~ TAXI" (Ref. 1)*has been’ 
collected for six commercial jet aircraft: These are‘the Boeing 707-320C, 727-200, 
and 747 and the McDonnell Douglas DC-8-63, DC-9-40 and the DC-10-10. The data °: 
presented in the form required by TAXI". A simulation was made for each aircraft 
taking off from two separate airfields, Dulles International in Washington D.C. 
and Will Rogers International Airport at Oklahoma City, Oklahoma. The purpose of 
using two profiles. being to point out the differences in aircraft response to 
different runway profiles. The calcomp plotted results are presented in the 
results section of this report. Geen: Gd Palins 
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This report alsv serves as a program users manual. /A sample problem simulating 
a Boeing 707-32%C during a constant speed taxi,(100 feet per second)> over the 
profile at the Will Rogers International Airport is included in the report. All 
of the fortran symbols used in "TAXI" are defined, and a complete listing of the 
program is~contained in appendices A, B, C, and D.’. Appendix E contains ;airplane 
data representing a “Typical” wide body tri-jet transport and two simulations 
using this data. 


ee 


~~ 
™ 


18, Distribution Statement 
This report is available through the 
National Technical Information Service, 
Springfield, VA. 22161, 


17, Key Words 
Runway Roughness 
Roughness 

Aircraft Response 
Runway Irregularities 


Airfield Pavement 
19, 


: 
+ ene 


21. No, of Feges 272, Puce 


Security Classif, ‘of this report) . Security Classif. (of thts page) 


Unclassifted Unclassified 137 
Form poT F 1700. i (8 n Reproduction of completed page outhorired iil 


iC 2 7 


Section 
I 
J! 


A 
1 
<= 


VI 
REFERENCES 


CEP ne 


APPENDIX A 


APPENDIX B 


APPENDIX C 


APPENDIX D 
APPENDIX E 


TABLE OF CONTENTS 


Page 
INTRODUCTION 1 
MATHEMATICAL MODEL c| 
EQUATIONS OF MOTION AND SOLUTION TECHNIQUE 10 
COMPUTER PROGRAM 14 
RESULTS 32 
CONCLUSIONS 76 

77 

COMPUTER LISTING FOR THE GENERAL COMPUTER 
PROGRAM CALLED "TAXI" 78 
COMPUTER LISTING OF "TAXI" USED TO SIMULATE 
CLASS E AIRCRAFT 93 


COMPUTER LISTING OF SUBROUTINE TAYLOR USED TO 
INCLUDE BEARING FRICTION IN THE STRUT FORCE 109 


FORTRAN SYMBOL DEFINITIONS 112 


TYPICAL WIDE BODIED TRI-JET TRANSPORT DATA 
AND SIMULATIONS 123 


Vv 


va ak AR ll ES EN 


i gated npr Se yt wb AP, “FP gda open zane awe 


nt Betty 


ie EE) 


std MY Gy He wt 


aN A 


LIST OF FIGURES 


Mathematical Model used in Taxi Simulation 

Typical Single Acting Oleo Pneumatic Landing Gear Strut 
Friction Force due to Wheel Drag 

Friction Force due to Nonperpendicular Orientation 
Friction Force due to Axle Offset 

Calcomp Plotted Results of Sample Problem 

Source Deck Setup for All Aircraft 

Class A Aircraft General Configuration 

Class A Aircraft Nose Landing Gear Tire Deflection Curves 
Class A Aircraft Main Landing Gear Tire Deflection Curves 
Class E Aircraft General Configuration 

Class E Aircraft Nose Landing Gear Tire Deflection Curves 
Class E Aircraft Main Landing Gear Tire Deflection Curves 
Class D Aircraft General Configuration 

Class D Aircraft Nose Landing Gear Tire Deflection Curves 
Class D Aircraft Main Landing Gear Tire Deflection Curves 
Class C Aircraft General Configuration 


Class C Aircraft Nose Landina Gear Tire Deflection Curves 


. Class C Aircraft Nain Landing Gear Tire Deflection Curves 


Class B Aircraft General Configuration 
Class B Aircraft Nose Landing Gear Tire Deflection Curves 


Class B Aircraft Main Landing Gear Tire Deflection Curves 


vi 


hi 


eines 


! 


a 


if 


bay MUU 


i ili A hl 


mM 


ul 


fi il 


CCR ecu 


— 


ii 


mh 


Figure 


23 


24 


25 


26 


27 


28 


29 


30 


3] 


32 


33 


34 


35 


36 


Plotted Results of a Class 
the Will Rogers Profile 


Plotted Results of a Class 
the Dulles Profile 


Plotted Results of a Class 
the Will Rogers Profile 


Plotted Results of a Class 
the Gulles Profile 


Plotted Results of a Class 
the Will Rogers Profile 


Plotted Results of a Class 
the Dulles Profile 


Plotted Results of a Class 
the Will Rogers Profile 


Plotted Results of a Class 
the Dulles Profile 


Plotted Results of a Class 
the Wil] Rogers Profile 


Plotted Results of a Class 
the Dulles Profiie 


Plotted Results of a Class B Aircraft During Taxi over 


A Aircraft Taking Off fron: 
A Aircraft Taking Off from 
E Aircraft Taking Off from 
E Aircraft Taking Off from 
D Aircraft Taking Off from 
D Aircraft Taking Off from 
C Aircraft Taking Off from 
C Aircraft Taking Off from 
B Aircraft Taking Off from 


B Aircraft Taking Off from 


the Will Rogers Profile with Strut Friction Included 


Plotted Results of a Class B Aircraft During Taxi over 
the Will Rogers Profile without Strut Friction Included 


Plotted Results of a "Typical" Wide Body Tri-Jet 
Transport Taking Off from the Will Rogers Profile 


Plotted Results of a "Typical" Wide Body Tri-Jet 
Transport Taking Off from the Dulles Profile 


Sainte 


: 


ise 


= 


a 


2 
= 
aS 
fy 
= 
me 
= 
Se 


ea 


o 


LIST OF TABLES 


aii ucla sae rena Oa i I 
y vO ner tet pet . 7 eae eaaletar eeeak ett) bi 
—! 


’ Table Page 
Aircraft Data for Conventional Aircraft 15 
Il Input Data Changes for Class E Aircraft 18 
Hl Runway Profile Magnetic Tape 19 
IV Profile Data of Will Rogers International Airport's a 
Taxiway E used in the Sample Problem 21 
a V Airplane Data for Sample Problem 29 
VI Typical Page of Listed Output Taken from the Results 
of the Sample Problem 31 
VII List of Input Data Fortran Names 35 
7 VIII List of Data Used to Simulate Class A Aircraft 36 
IX List of Input Data Fortran Names 41 
X List of Data Used to Simuiate Class £ Aircraft 42 
XI List of Input Data Fortran Names 47 
XII List of Data Used to Simulate Class D Aircraft 48 
a XIII List of Input Data Fortran Names 53 
4 XIV List of Data Usec to Simulate Class C Aircraft 54 
i: XV List of Input Data Fortran Names 59 
¥ XVI List of Data Used to Simulate Class B Aircraft 60 
XVII List of Data Used to Simulate the Wide-Body 


Tri-Jet Composite 126 


vili 


LIST OF SYNSOLS 


Definition Unit 
Distance from rear main landing gear to aircraft CG a 
Coefficients calculated to fit a polynomial curve to NL 
a runway Segment 
Hydraulic piston area in? 
a Ay Effective orifice area int 
Ht B Distance from nose landing gear to aircraft CG in 
ED Cc Distance from front main landing gear to aircraft CC in 
c Damping coef fifient 1/sec z 
: C, Orifice coefficient ND q 
Fa Pneumatic strut force Ibs = 
Fan Total airplane aerodynamic drag Ibs 
Fy oFy Landing Gear Drag and Normal Forces : Ibs 
Fe Friction force Ibs 
Fee Strut seal friction ibs 
’ Fr, Hydraulic strut force Ibs Z 
| FopefsaFs3 Total landing gear strut forces Ibs 4 
4 Fe oF ey FoF eg Landing gear tire forces Tbs q 
4 Fy Total airplane thrust Ibs E 
4 . Frp Main landing gear tire drag force Tbs J 
i 7 FoFy Strut bearing friction forces, upper and lower Ibs z 
q g Gravitational! acceleration constant in/sec™ 
q ne Total airplane pitching inertia Ibs/in secé 
q k Linear tire spring constant Ibs/in 


| 
ST 


ix 


at 


i 


NA 


iit 


ta ad 


Definition 

Total airplane lift 

“iigt27g° "3/9 

Airplane mass 

The generalized mass for the ith mode 
Fully extended strut air pressure 


Def ‘ection due to bending 


rray of nondimensional time dependent coordinates 
whie ) weigh the amount of motion due to the ith mode 


in the total motion of the aircraft 
Time devivative of q: 
Time derivative of q; 

Strut displacement 

Time derivative of § 

Second time derivative of § 

Taylor series sample solution time step size 
Tire deflection 

Fully extended strut volume 

Total airplane weight 

Unsprung landisg gear weights 

Coordinate for sample Taylor series solution 
Time derivative of x 


Time derivative of x 


Horizontal translation of the aircraft down the runway 


Time derivative of X 
Time derivative of X 


Vertical displacement of the vehicle CG 


Time derivative of Z 


in/sec 
in/sec” 


sec 


in/sec 
in/sec* 
ft 
ft/sec 
ft/sec* 


in 


in/sec 


Definition 
Time derivative of Z 


Vertical displacement of the unsprung mass of gear 
number |, 2 and 3 


Time derivative of Z) 129 523 


Time dsrivative of 2 sty25 


histance from main gear axle to the waterline of 
{%: aircraft CG 


Mocal deflection of the ith mode at gear 1, 2 and 3 
Aircraft pitch angle 

Time derivative of 6 

Time derivative of 6 

Hydraulic fluid density 

Percent of critica! damping factor for the ith mode 


Sliding coefficients of friction for strut lower 
and upper bearings 


Natural frequency of the ith mode 


Unit 
é 2 
in/sec 


in 


in/sec 


in/sec” 


rad/sec 


rad/secZ 


sligs/in? 
ND 


ND 


rad/sec 


ony 


ng " i i i 


ES REDLEAF ORER SAD THE WOH Pg EO 


PUEBLA 


FE 


FEE a aR 


uu, UNE) GU 


a 


4 
43 
= 
EB 


SECTION i 


INTRODUCTION 


Very early aircraft operated from grass airfields or gr 
untreated runways. However, with the development of more sophis® 
aircraft requiring higher and higher ground speeds, runway surfa¢ 
became of increasing concern. Today’ s modern aircraft require a pave- 
ment system structurally capable of supporting gross loads approaching 
a million pounds. In addition, the pavement surface musi provide adequate 
skid resistance and must be of sufficient smoothness to prevent undue loss 
of fatigue life in the aircraft. The response of aircraft to the pavement 
irregularities has become of increasing concern in recent years. A great 
deal of emphasis has been placed upon the comfort and safety of aircraft 
ground operations and other factors related to the interaction between 
the pavement and the aircraft. 


Because the runway and the aircraft form a coupled system with the 
pavement profile providing a displacement input that can radically affect 
the behavior of the vehicle, there is an urgent need for a means to measure 
pavement roughness. Yhe Federal Aviation Administration has tong been 
aware of this need, and in 1973 initiated research and development efforts 
to quantify pavement roughness. This effort began with the study of a 
means for rapidiv collecting the runway profile. As a result of that 
effort, a profiling system using a laser beam as a horizontal reference 
datum was developed (Reference 1). 


The next step was to establish roughness criteria. Considerable 
effort had been expended in this area, most of it dealing with the sub- 
jective rating of ride quality. However, ride quality did not deal directly 
with the pavement, which was the factor causing unacceptable aircraft 
response. Therefore, in its efforts to define runway roughness, a statis- 
tical approach to the runway profile itself was taken (Reference 2). This 
approach dealt with the root means square of the deviation of the profile 
about a normalized line. From this comes criteria for evaluating a runway 
profile as being acceptable as long as the RMS level is below 0.32 inches. 
Marginal conditions are equated to RMS levels of 0.32 to 0.36 inches and 
any runway having a RMS level of greater than 0.36 inches is considered 
to be unduly rough. 


As a check of this criteria and as a means of evaluating the simu- 
lated repair of a runway having undesirable roughness, an existing taxi 
code developed for the US Air Ferce was modified to provide a way for 
calculating the dynamic response of commercial aircraft to runway irre- 
yularities. The modification of that computer code is the subject of 
this report. The use of this computer program is not to compare ride 
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quality between various classes of aircraft but rather to be used by 
the operating airlines and airport owners and operators as a means 
of identifying rough pavement areas..and to help insure the timely 
maintenance of pavements in crder to prevent the develcpment of undue 


roughness. 


Details of the development of the computer code modified in this 
report is provided in Reference 3. The development of the computer 
code is only briefly summarized in this report. in order te make the 
computer code useful to al] commercial aircraft. these aivcraft were 
divided into five classes based on gross weight as shown below: 


Aircraft Classes 


Class A Aircraft Less than 150,000 pcunds 
Class B Aircraft 158,600 - 300,000 pounds 
Class € Aircraft 300,000 - 450,000 pounds 
Class 0 Aircraft 450,000 - 600,000 pounds 
Class E Aircraft Greater than 600,000 pounds 
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SECTION Ii 
MATHEMATICAL MODEL 


The general airplane/runway mathematica] medel used for this study 
was the basic mathematical model developed in Reference 1. A detailed 
description of the components that make up this general model, as well as 
the assumptions made are shown in Reference 1. This report presents, in 
summary form, the landing gear strut and tire representation, the airplane 
rigid body and flexible body representation, the runway profile represen- 
tation, the equations of motion and the solution technique. 


The general model is represented as a symmetrica! body with either a 
nose and single strut main landing gear typical of class A through D air- 
craft. Each landing gear strut is assumed to have point contact with the 
profile and it is assumed that the main landing gear traverses the same 
profile as the nose landing gear except at a later time. The model has 


aerodynamic lift and drag and thrust applied at the aircraft's center of 
gravity. 


The fuselage is free tu pitch, plunge and translate horizontally down 
the runway. In addition to these rigid body degrees of freedom, up to 15 
flexible modes of vibration are included. This fuselage motion is con- 
trolled by the landing gear strut forces, lift, drag, thrust, and the 
resisting parameters of aircraft mass and inertia. Figure 1 is a schematic 
drawing of the mathematical model used for this simulation. 


The landing gear struts are non-linear, single acting oleo-pneumatic 
energy absorbing devices (Figure 2) and are represented in the model as 
the sum of the three forces, penumatic, hydraulic, and strut bearing 
friction forces. The pneumatic force, which is the laraest of the three 
is represented by the equation: 


Fras (1) 


where: 


pneumatic strut force 


the fully extended strut pressure 


the fully extended strut volume 


Nt 


the pneumatic piston area 
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strut stroke 
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The hydraulic or damping force is given by the equation: 
PAR? SISI (2) 
N24 Ae 
where: 
Fs = hydraulic strut force 
P, = density of the hydraulic fluid 


A, = the hydraulic piston area 


A_ = effective orifice area (constant orifice minus 
metering pin area) 


Cy = orifice coefficient (use 9.9) 


fre 
" 


strut piston velocity 


The third strut force is the strut bearing friction force and is 
included in the model only when articulated or asymmetrically loaded 
struts are being simulated. In symmetrically loaded struts the friction 
force is neglected (si: Reference 1). The following derivation of the 
friction force was taken from Reference 2. 


The force on upper and lower strut bearings due to rolling friction is 
shown in Figure 3: 


*D, = .025 Fy 
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Figure 3. Friction Force Due to Wheel Drag 
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where: 
Fo = landing gear drag force 
Fy = landing gear normal force 
"FF, Fr = bearing forces at the upper and lower bearings 


7 respectively required to balance the lateral 
loading on the piston 


oo) 
u 


the piston length 


a 
" 


b the strut fully extended bearing separation 
s = the strut stroke 


The force on upper and lower strut bearings due to nonperpendicular 
orientation is represented by Figure 4: 
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Figure 4. Friction Force Due to Nonperpendicular Orientation 
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Finally, as represented in Figure 5, the force on the upper and ltewer strut 
bearings due to axle offset is: 


Figure 5. Friction Force Due to Axle Offset 


For a landing year strut having wheel drag, an angle from the vertical 
and an axle offset, the total friction force would be as follows: 


Fe = tu, Fe +u,, Fh + Fel Ist 


ad 


$s 
where: 
Fe = friction force 
My? bp F coefficients of sliding friction at the upper and 
. lower bearings, respectively; u, = 0.10, y, = 0.15 
(“he lower bearing is generally not as well lubricated 
as the upper bearing. ) 
Fe = "strut seal friction at zero lateral loading 
0 
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aa a 
H The tire force is represented by the linear equation 
Fe = k Ty (3) 
where: 
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t tire force 


onl 
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D tire deflection 


age 


a” 
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linear tire spring constant 


The runway elevation data is input into the model in two-foot 
increments, three data points at a time for each landing gear as required. 
4 The profile is made continuous by fitting the following polynomial throush 


Hut 


the three elevation data points and the slope at the end of the previous 
profile segment: 


V(x) = ay + ax + ax" + a,x° (4) 
where: 


21 2,3,4 * coefficients derived from the elevation and slope data 
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SECTION III 


rip 


Equations of Motion: 


notation shown in Figure 1. 


Rigid Body Equations of Motion 


4 2 = (Foy + Foo + Fe3 #+tL- W)/Meg [c.g., vertical acceleration] (5) 
Zy= (Fey > Fey > Wy) /M) (6) 
- [unsprung landing gear j 
29= (Fro ee Wo) /M, mass acceleration} (7) | 
ee ae al ae - 
Gx (Fy A+ Feo B+ Fro’) - Fo3 Os {pitching acceleration] (9) 
= = = = f > * . 
Xx (Fy Feo Fap/ clam) [horizontal translation acceleration] (10) 
where: 
4 Foy. Foo» F53 = total landing gear strut forces 


Fire Feo Fig = tire forces 


= 
r) 


a Wis Wos W3 = unsprung landing gear weights 


A, B, Ey» = feet arms 
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EQUATIONS OF MOTION AND SOLUTION TECHNIQUE 


The differential equations of motion for the mathematical model were 
derived by application of Hamilton's principle of minimum energy. The 
general form of these equations is shown below and corresponds to the 


aircraft mass, weight and pitching inertia 


Frp? FD = }if>. thrust and tire and aerodynamic drag forces 


d Fap act through the center of gravity] 
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Flexibility Equations of Motion 


- : : 
May “Ei1Fs) + Si2Feo + §i3F53 - 25;4)M;q; - 4 M;q; for the ith mode (11) = 
where: | 
M. = the generalized mass F 
Giy» Sig» 63 = modal deflections at gear location 1, 2 and 3 i 


w, = modal frequency 


i 


&; = damping factor 
4° Ws a = generalized coordinate and time derivatives 


The sign convention corresponding to Figure 1 and the equations 
of motion is as follows: 


Z = Vertical Displacement + up 
a 6 = Pitch + nose down 
7 q = Deflection Due to Bending + up 
: X = Horizontal Translation + forward 


Solution Algorithm: 


The technique used for solving the coupled nonlinear differential 
equations of motion that describe the simulated aircraft is a three- 
term Taylor series. For example, the equation: 


R = -cx -kx 


i (72) 

4 The three-term Taylor series representation can be written as: 

s x x (at)? 13 
X ret) 7 % + * ay (at) + x (1) —— (3) 


ca 


where: I = | oN 


The values for x, x and x from the previous step are substituted inio 


Lh ACES 
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equation (13) and a new value for x is obtained. Differentiating 
equation (13) we obtain for the velocity x, the expression 


Keren) = ry + Ryle) (14) 


The values for x and x are then substituted into equation (14) 


and a new value of x is found. 


This entire process is repeated with 


the new values of x and x to obtain the next point in the solution. 


In Section IV Table Il identifies and shows the nodifications of 


cards which are necessary for a Class E aircraft. 
remain the same. 


Table III shows the format for inputing the runway 


profile data. 


The sample problem at the end of this section simulates a Class 


C€ aircraft traversing the pr 
International Airport at a constant speed of 100 feet per second. 


ofile cf Taxiway E at the Nill Rogers 


sample problem output in order of occurrence is as follows: 


Runway Profile Listing 

List of Input Data Used 

One Page of Listed Time History Output 
Time History Plot 
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All other data cards 
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SECTION IV 


COMPUTER PROGRAM 


The computer preqram, used to predict the dynaniic response of the 
commercial aircraft described in the text section, is the basic program 
described in Reference 1 with a few minor changes. These program changes 
include the capability to simulate strut hearing friction forces, multiple 
main landing ‘struts, root mean square calculations, compressed plot 
capability. co 


Appendix A contains a complete listing of the general computer program 
for all aircraft except the Class — aircraft. A complete listing for the 
program used to simulate the Class E aircraft is shown in Appendix 8B. To 
Simulate a Class B aircraft with tne inclusion of strut bearina friction, 
the separate “TAYLOR® subroutine listed in Appendix ©, should be substi- 
tuted for the "TAYLOR" subroutine included in the general program. If it 
is desired to run Class B aircraft without strut bearing friction, the 
general program can be used. Appendix 5 contains a complete listing of 
the FORTRAN symbols used. 


The following tables I through V° contain the form of the aircr.ft 
input data required for an aircraft airfield simulation. The data cards 
are sequenced as they must appear in whe aircraft data deck. Fur each data 
card, variable names, definitions, un-ts. card columns, and format field 
Specifications are given. 


Table I contains the form of tre aircraft data deck for all con- 
ventional (Class A through D) landing gear aircraft. These aircraft all 
have main landing gear struts which a:¢ single acting and are non-articulated. 
Table I also contains the data used in the sample problem at the end of 
this section. 


ressec Or suppressed j rinting, and other minor improvements. 
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Card 
Column 


Card | 


1-80 


Card. 2 


4-10 

11-20 
21-30- 

31-42 


Card 3 


1-10 
11-20 


Card 4 


1-10 
11-20 
21-30 


Card § 


1-10 
11-20 
21-30 


Card 6 


1-10 
11-20 
21-30 
31-40 


Format 


8A10 


FI0.1 
F10.1 
F1O.1 
FI2.0 


FIO.2 
FI0.2 


Fl0.2 
FI0.2 
FI0.2 


FIO.4 
FIC.4 
FIO.4 


FI0.2 
F10.2 
FI0.2 
F10.2 


TABLE | 


AIRCRAFT DATA FOR CONVENTIONAL AIRCRAFT 


Section | (cards 1-5) - General Airptane Data 


Variable 
Name 


PLANE 


W 
A 
8 
MMI 


PSARM 
TAILRM 


SPEED 
THRUST 
TAKEOFF 


CL 
AREA 
cD 


Data for 


Sample 


Problem 


Class C 


678. 


84700000 


757 
637 


100 
48000. 
289. 


-603 
2890. 
-03 


Definition 


Airplane Being Simulated 


Vehicle Weight (1b) 

Distance Main Gear to CG (in) 
Distance Nose Gear to CG (in) 

Mass Moment of Inertia (Ib in sec*) 


Distance of Pilot Station to CG (in) 
Distance of Tail Station zo CG (in) 


tnitial Velocity of Airplane (ft/sec) 
Total Airplane Thrust (1b) 
Airplane Ratation Speed (ft/sec) 


Lift Coefficient 
Wing Area (ft2) 
Drag Coefficient 


Section 2 (cards 6-11) -_ Main and Nose Gear 


1659. 
432. 
2. 

'. 


Unsprung Weight of Each Main Gear (1b 
Unsprung Weight of Nose Gear (Ib) 
Number of Main Gear Struts 

Number of Nose Gear Struts 


aL NTT 


se 


fi : igs 04 


Weg wae EAL FB 


POLO NE 


SY 


Card 


1-10 
11-20 
21-36 
31-40 


Card 8 
i-10 
11-20 
21-30 
31-40 
4\=50 
51-60 
Cars 9 
1-10 


11-20 


a 


i-10 
11-20 


1 


1-10 


7 Card 13 


Column 


Card 10 


Card |1 


Card i2 


FI0.5 


15 
15 


Variable 


Name 


TSN 


DX 


1FPLOT 
(FLIST 


Data for 
Sample 
Problem 


oT 


13.91 
19.64 
66.80 
78.47 


265. 
243. 
335. 

1816.6 
3.14 
1.23 


91.0 
91.0 


25050. 
13000. 


-00) 


Definition 


Hydraulic Pisten Area Nose (in*) 
Pneumatic Piston Area Nose fin?) 
Hydraulic Piston Area Main (in?) 
Pneumatic Piston Area Main {in@} 


Nose Strut Préioad Pressure (1b/in?) 

Main Strut Preload Pressure (Ib/in2) 
Fully Extended Nose Strut Air Volume (in?) 
Fully Extended Main Strut Air Volume (in3) 
Orifice Area Main (in?) 

Orifice Area Nose (in*) 


Distance from Axle to C& Waterline 
Main Gear Strut Unloaded (in) 
Distance from Axle to CG Waterline 
Nose Gear Strut Unloaded (in) 


Main Tire Spring Constant Per Strut (1b/in 
Nose Tire Spring Constant Per Strut (1b/in 


Integration Step Size 


O-Plot; I-No Plot 
C-List; I-No Plot 


Section 3 (cards 13-16)-Metering Pin Description 


Number of Slope Changes Nose Gear 


Data for 


Card Variable Sample 
Column Format Name Problem Definition 


f 


*Card 4A, 14B,.... 


1-10 F10.3 STROKN (1) * Stroke Corresponding to Metering Pin 
Diameter, Nose Gear 


11-20 F10.3 PINDN (1) * Metering Pin Diameter, Nose Gear (in) 
Card 15 


1-5 15 NSCM :) Number of Slope Changes Main Gear 
| 


*Card 16A, 16B,.... 


1-10 F10.3 STROKM (1) * Stroke Corresponding to Metering Pin 
Diameter, Nose Gear 


11-20 F10.3 PINDM (1) * Metering Pin Diameter, Main Gear (in) 


Section 4 (cards 17-19)-Flexibility Data 
Card 17 


1-5 15 NEM 10 Number of Flexible Modes 


**Card 18A, 18B..... 


1-10 FI0.3 SIMAIN (1) we Mode Shape Deflection Main Gear 
11-20 FI0.3 SINOSE (1) ax Mode Shape Jeflection Nose Gear 
21-30 F10.3 SicG (%) ae Mode Shape Deflection CG 
31-40 F10.3 SITAIL (1) kr Mode Shape Deflection Tail Station 
41-50 F10.3 sips. (1) wet Mode Shape Deflection Pilot Station 
**Card I9A, 19B,.... 

1-15 FIS.2 gy (1) we Generalized Mass (Ib sec*/in) 
16-25 F10.3 OMEGA (1) wx Modal Frequency ({rad/sec) 


* 


*0ne card ¢s required for each stroke-metering pin combination read into the 
program. 


il 


Data for Sample Problem is shown in Table V. 


é 


**0ne card is required for each flexible mode. Data for Sample Problem is shown 
in Table V. 


vpenpe sg 7 7 1 bie ian ne 
sil Hata etn abhi A a i Lavage aa lg aig hil i 


MW 


Card 
Column 


Card 2 


1-10 
11-20 
21-390 
31-40 
4} -52 


a ss Pc cL 


Card 6 


1-10 
11-20 
21-30 
31-40 
41-50 


Card 7 


1-10 
11-20 
21-30 
31-40 
4}-50 
51-60 


AU 


Card 8 


uly 
4 


HHO 


1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 
71-80 


nb 


HAM aaa 


7 Hy 
tt 


Card 10 


1-10 
11-20 
21-30 
31-40 


Format 


FIO. 
F10. 
FIO. 
FIO. 
Fl2. 


Cow me ee ee 


FIO. 
FIO. 
F10. 
FIO. 
FI0,2 


BOR RN PO 


F10.5 
F10.5 
F10.5 
FIO.5 
FI0.5 
FI0.5 


FI0.5 
FIO.§ 
FI0.5 
FI0.5 
FI0.5 
FIO.5 
FIO.5 
FIO.5 


FIO. 
FIO.1 
FI0.1 
FIO.1 


TABLE 11 


INPUT DATA CHANGES FOR CLASS E AIRCRAFT 


Variable 


Name 


Definition 


Vehicle Weight (Ibs) 

Distance Front Main Gear to CG (in) 
Distance Nose Gear to CG (in) 
Distance Rear Main Gear to CG (in) 
Mass Moment of Inertia (Ibin sec‘) 


Unsprung Weight of Each Front Main Gear (1bs) 
Unsprung Weight of Nose Gear (Ibs) 

Number of Main Gear Struts 

Number of Nose Near Struts 

Unsprung Weight of Each Rear Main Gear (Ibs) 


Hydraulic Piston Area Nose (in?) 
Pneumatic Piston Area Nose (in*) 
Hydraulic Piston Area Front Main (in@) 
Pneumatic Piston Area Front Main (in“) 
Hydraulic Piston Area Rear Main (in‘*) 
Pneumatic Piston Area Rear Main (in) 


Nose Strut Preload Pressure (Ibs/in?) 

Front Main Strut Preload Pressure (Ibs/in2) 
Fully Extended Nose Strut Air Volume (in3) 

Fully Extended Front Main Sgrut Air Volume (in?) 
Orifice Area Front Main (in*) 

Orifice Area Nose (in2) 

Rear Main Strut Preload Pressure (Ibs/in?) 

Fully Extended Rear ‘Jain Strut Air Volume (ind) 


Front Main Tire Spring Constant Per Strut (ibs/in) 
Nose Tire Spring Constant Per Strut (Ibs/in) 

Rear Nain Tire Spring Constant Per Strut (ibs/in) 
Orifice Area Rear Main (in2) 
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Runway Profile Magnetic Tape : 


The runway profile is read into the program from a magnetic tape or 
permanent file. The format for this tape is shown in Table III. 


TABLE Jt! 


RUNWAY PROFILE MAGNETIC TAPE 


Variable 

Column Format Name Definition 

Read 1 

1-80 8A10 SITE Runway Profile and Direction 
Read 2 

1-6 16 NPTSS Number of Runway Elevation Points 

" #Read 3, 4... .N+2 
1-70 10F7.3 ELEV Runwey Profile Data 


*One read required for every ten runway profile elevation points. 


Figure 7 contains a schematic diagram of the source deck setup for 
all aircraft simulations. 


In order to simulate Class B aircraft with the strut friction force 
included, changes must be made to the basic program. This change is 
required as a result of the main and nose landing gear having a sianifi- 
cant angle from the vertical and an axle offset on the main landing gear. 
The Class B aircraft source deck is formed by removing the subroutine 


Taylor from the basic source deck and replacing it with the Taylor shown 
in Appendix C. 
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- Source Deck Setup for All Aircraft 
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SECTION V 
RESULTS 
The mathematical model and computer program presented in this report 


is capable of simulating five classes of commercial jet aircraft. Classi- 
fication is made on the basis of gross weight as follows: 


Class A Less than 150,00C 1b 
Class B 150,600 - 300,000 Ib 
Class C 300,000 - 450,000 1b 
Class D 450,000 - 600,000 1b 
Class E Greater than 600,000 1b 


The data used to represent these aircraft and typical simulations using 
the data for each class of aircraft is presented in this section. 


The data set for each aircraft class contains a general configuration, 
(Figures 8, 11, 14, 17, and 20), a description of input variable fortran 
names, (Tables VII, IX, XI, XIII, and XV), two pages of input data (Tables 
VIII, X, XII, XIV, and XVI), and a set of tire deflection curves, (Figures 
9, 10, 12, 13, 15, 16, 18, 19, 21, and 22), for both main and nose landing 
gears. The first page of input data contains al? of the rigid body air- 
plane and landing gear data. The second page of data contains the aircraft's 
modal data and the aircraft's initial conditions as computed by TAXI. 


The aircraft data presented was used to simulate each aircraft taking 
off from two separate profiles, Will Rogers International and Dulles 
International Airports. The calcomp plotted results are shown in Figures 
23 through 32. It is evident that all airplanes responded to the prominent 
bump which occurred at approximately 2700 feet on the Will Rogers profile, 
although some to larger degree than others. With the exception of this 
one bump it is difficult to correlate runway roughness to the response 
of all aircraft. Some airplanes responded to a certain section of a 
profile while others did not. Some airplanes responded predominately in 
pitch such as the Class C and D aircraft while others responded in pitch 
and plunge such as the Class A and B aircraft. This points out the signi- 
ficance of such parameters as gear stiffness, degree of coupling between 
the rigid body i:-.es, gear spacing and others. 


For ali aircraft classes except Class B the rms value of the Pilot's 


Station vertical acceleration was less over the Dulles profile than that 
of the Will Rogers profile. For all aircraft classes except Class B and 
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D, the center of gravity vertica} acceleration rms level is less for the 
Dulles profile. For these three exceptions, the rms level was the same 
for both profiles. From a purely response standpoint then, it can be 
said that the Dulles International runway profile is smoother than the 
Will Rogers International orofiie. 


Figures 33 and 34 show the difference in response of Class B air- 
craft with ana without strut bearing friction included in the simulation. 
The simulation for both runs were at a constant speed of 100 feet per 
second over the Will Rogers profile. The center of gravity rms level 
was increased with the inclusion of strut friction while it was reduced 
slightly at the pilot's station. The jerky motion evidenced by the 
response at the cg was expected and was caused by the constant "locking 
up" of the struts due to the binding forces on the upper and lower bear- 
ing surfaces of the landing gear struts. It is more evident at the cg 
because most of the aircraft's weight, and resulting strut forces, is 
on the main landing gear which is closest to the aircraft's cg. 
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Figure 11 - Class E Aircraft General Configuration 
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SECTION VI 
CONCLUSI “NS 


Data have been collected and reduced tc the form required by the 
computer program "TAXI" for simulation of five classes of commercial jet 
aircraft as follows: 


Aircraft Class Gross Weight Range-Pound 
Class A Less than 150,000 

Class B 150,000 - 300,000 

Class C 300,000 - 450,000 

Class D 450,000 - 600,000 

Class E Greater than 600,000 


Each of these aircraft has been successfully simulated traversing 
the profiles of Will Rogers International Airport in Oklahoma City, 
Oklahoma and Dulles International Airport in Washinoton DC. Results of 
these simulations indicate that aircraft response to runway roughness 
is highly dependent on aircraft parameters such as gear stiffness, gear 
spacing and the degree of coupling between the rigid body modes of 
vibration, as wel] as the cegree of runway roughness. 


The addition of the strut bearing friction forces in the simulation 
of Class B aircraft has a small effect on the vertical acceleration level 
at the pilot's station. 
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APPENDIX A 


COMPUTER LISTING FOR THE GENERAL COMPUTER PROGRAM CALLED "TAXI" 
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APPENDIX C 


COMPUTER LISTING GF SUBROUTINE TAYLOR USED TG INCLUDE 
BEARING FRICTION IN TRE STRUT FORCE 
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APPENDIX D 


FORTRAN S“MBOL DEFINITIONS 


This appendix contains an alphabetical listing of the Fortran 
variables used in the program "TAXI" categorized by the subroutine in 
which they are defined. In cases where a variable is defined in two 
or mere Subroutines, it is listed under the subroutine in which it is 
used most often. Sceme symbols used in the Class E Aircraft computer 
code are not listed. These variables are those which have been formed 
by adding an R to a variable name which is contained in the basic 
"TAXI" computer code. The R refers to the rear set of main gear. 


FORTRAN SYMBOL DEFINITIONS 


3 TAXI 
i co 
H SYMBOL DEFINITION 
5 
A Distance from CG to front main gear 
AAM Pneumatic area, main gear 
AAN Pneumatic area, nose gear re 
AHM Hydraulic area, main gear ie 
AHN Hydraulic area, nose gear ey 
AM Coefficient of polynomial fit to runway profile segment, iS 
rear main gear i 
AN Coefficient of polynomial fit to runway profile segment, 
nose gear 
AREA Aircraft wing area 
B Distance from CG to nose gear 
BM Coefficient of polynomial fit to runway profile segment, 
rear main gear 
BN Coefficient of polynomial fit to runway profile segment, i 
nose gear E 
4 C Distance from CG to rear main gear {Class E simulation only! Ps 
4 cD Coefficient of drag d 
CGACC Array containing CG accelerations for Calcomp plot 3 
CGOUT Total CS acceleration 3 
cL Coefficient of lift 3 
q cM Coefficient of polynomial Fit to runway profile segment, 
a rear main gear 
CN Coefficient of polynomial fit to runway profile segment, 


nose gear 
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DODPLOT 


DiSTAN 


DM 


DN 


DX 


02 


ELEV 
ELEV! 
ENDRUN 


GM 


HDR 


E 1FPLOT 


IFLIST 


I VAL 


IRVAL 


me 


« 
mh 


IXLONG 
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SYMBOL 


DEFINITION 


Array containing a.ccraft distance down the runway for 
Calcomp plot 


Distance down the runway used in removing overall slope 
from runway profile 


Coefficient of polynomial fit to runway profile segment, 
rear main gear 


Coefficient of polynomial fit to runway profile segment, 
nose gear 


Time step for integration 


Variable which comeendates for the overlap of two adjacent 
runway segments 


Array containing runway profile elevations 
Elevation of first runway profile point 
Length of runway 


Array containing generalized masses for each flexible 
mode of viuration 


Counter for printing header on printed output 
Index variable 


Var‘able which contains decision to produce plotted output 
or not 


Variable which contains decision to produce a listed output 
or not 


Subscript variable for runway markers on Calcomp plot 


Integer truncation of distance between nose and rear main 
gear 


integer truncation of distance between rear main gear and 
front main gear (Class E simulation only) 


Integer truncation of length of time axis on Calcomp plot 


Subscript variable for runway profiie 
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Symoo! 


LD! 


x 
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DEFINITION 
Counting variable for runway profile input 
Counting variable for runway profile output listing 


Subscript variable for storage of pilet station 
acceleration time history 


Subscript variable for storage of aircraft speed and 
distance for Catcomp plot 


Runway distance for runway profile listing 

Counting variable for runway profile input 

Counting variable for runway profile listing 

Counting variable for printing out output header first time 
Mass of entire aircraft 

Integer truncation of totai simulation time 

Mass of unsprung portion of one main landing gear 

Pitching moment of inertia about aircraft center of sravity 
Mass of the unsprung portion of the nese landing gear 
Length of runway divided by 1000 feet 

Number of flexible modes 

Subscript variable for CG acceieration time history 

Number of runway profile points 

Number of slope or area changes on main strut metering pin 
Number of slope or area changes on nose strut metering pin 
Defines run as taxi or takeoff 

Index variable for normalization of runway profile 

Area of orifice hole, main gear 


Area of orifice hole, nose gear 
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SYMBOL 
OMEGA 
PAQM 
PAON 


PINDM 


PINDN 


PLANE 


PP 


PROF 
PROF IO 


PSA 


QD 
QDD 
QDOCG 
QDOPS 
QODTAL 
REACTM 
REACTN 


RM 


DEFINITION 

Array of flexible mode frequencies 

Preload pressure of main gear strut 

Preload pressure of nose gear strut 

Array containing main gear metering pin diameters for con- 
ventional aircraft and net orifice areas for aircraft with 
metering tubes or flutted metering pins 

Array containing nose gear metering pin diameters for con- 
ventional aircraft and net orifice areas for aircraft with 
metering tubes or flutted metering pins 


Aircraft being simulated 


Counting variable for storage of distance and speed for 
Calcomp plot 


Rurway profile time history elevations 

PROF (NN+2)X10 

Pilot station acceleration 

Array containing pilot station acceleration time history 
Distance from pilot station to CG 

Array of non-dimensional time dependent coordinates which 
weight the amount of motion due to each flexible mode in 
the total motion of the aircraft 

Time derivative of Q 

Time derivative of QD 

CG acceleration due to flexible motion 

Pilot station acceleration due to flexible motion 

Tail station acceleration due to flexible motion 

Static, total force at main gear 


Static, total force at nose gear 


incremented variable for determining position of runwav 
markers 
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SYMBOL 
RMARK 
SicG 
SIMAIN 
SINOSE 
SIPS 
SITAIL 
SITE 


SLM 


SLN 


SLP 
SPEED 
SSPLOT 
STORE| 
STORE2 
STORE3 
STORES 
STROKH 


STROKN 


SXM 
SXN 
TAILAC 
TAILRM 


TAKOFF 


DEFINITION 
Array containing runway markers positions 
Mode shape deflection of CG 
Mode shape deflection at main landing gear 
Mode shape deflection at nose landing gear 
Mode shape deflection at pilot station 
Mode shape deflection at tail station 
Location of runway 


Distance from CL of main gear axle to CG of aircraft 
with strut fully extended 


Distance From CL of nose gear axle to CG of aircraft 
with strut fully extended 


Overall slope of runway profile 

Initial speed of aircraft 

Array of velocity of aircraft for Calcomp plot 
Temporary storage space for CG accelerations 


Temporery storage space for CG accelerations 


Temporary storage space for pilot station accelerations 
Temporary storage space for pilot station accelerations 


Array of strut siroke values corresponding to metering 


pin values (PINDM), main gear 


Array of strut stroke values corresponding to metering 


pin values (PINDN), nose gear 
Number of main gear struts 
Number of nose gear struts 
Acceleration at tail station 
Distance from tail station to C6 


Rotation velocity of aircraft 
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SYMBOL 
THRUST 


TIME 


TIMEX 


TIMEI 


TSM 
TSN 
TYPRUN 


vom 


XLONG 
XLONG2 
XPROF 
XSTOP 


YP 


YPR 


YPN 


ZRDOT 


Total thrust of aircraft 


Array of simulation times at which CG acceleration time 
history points are stored 


Counter variable for printed output 


Array of simulation times at which pilot station time 
history are stored 


Tire spring constant, main gear 

Tire spring constant, nose gear 
Defines simulation as takeoff or taxi 
Main gear strut fully extended volume 
Nose gear strut fully extended volume 
Weight of aircraft 

Main gear unsprung weight 

Nose gear unsprung weight 

Simulation time 

Length of tine axis for Calcomp plot 
XLONG/2 

Location for printing of NOSE GEAR TRACK" on Calcomp plot 
XLONG+5 


Array containing runway segment elevation points and slope 
from end of previous segment, rear main gear 


Array containing runway segment elevation points and slope 
from end of previous segment, rear main gear (Class E 
simulation only) 


Array containing runway segment elevation points and slope 
from end of previous segment, nose gear 


Slope of runway segment at end point, rear main gear 
(Class E simulation only) 


SYMBOL. 


ZNDOT 


TAYLOR 
AOM 
AON 

“BSN 


*BSM 


COM 
a CON 
DRAGA 


DRAGT 


*FFON 
“FFOM 
“FUN 


PLM 


FSM 
FSN 


4 FSTN 


FSTI 


owt 


FSTZ 


ye 
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3 FTN 
*FUN 
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*FUM 


DEFINITION 


Slope of runway segment at end point, nose gear 


Net orifice area, main gear (OAM~metering pin area) 
Net orifice area, nose gear (OAN-metering pin area) 
Nose landing gear minimum bearing separations 

Main landing gear minimum Dearing separations 
Damping coefficient, main gear 

Damping coefficient, nose gear 

Aerodynamic drag 

Rolling Drag 

Nose landing gear seal friction 

Main landing gear seal friction 

Nose ianding gear lower friction force 

Main landing gear tower friction force 

Total force in all main gear struts 

Total force in all nose gear struts 


Net force on secondary piston, nose gear (Class E 
simulation only) 


Net force on secondary piston, rear main gear 
(Class E simulation only) 


Net force on secondary piston, front main gear 
(Class E simujation only) 


Force in tires, main gear 
Force in tires, nose gear 


Nose landing gear upper friction force 


Main landing gear upper friction force 


*Class B simulation with Strut Friction only 


SYMBOL 
*OF FEN 
*OFFM 
*PHIN 
“PHIM 
“PLN 
*PLM 
QTDM 


QTON 


QTM 
QTN 


SLOPEM 


SLOPEN 


- SSM 


=E 
i 
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) TD{2) 
TO(4) 
TD (6) 
TD(8) 


DEFINITION 
Nose landing gear axle offset 
Main landing gear axle offset 
Nose landing gear angle from the vertical 
Main landing gear angle from the vertical 
Nose landing gear piston length 
Main landing gear piston length 
Total velocity due to flexible modes at main gear 
Total velocity due to flexible modes at nose gear 
Total deflection due to flexible modes at main gear 
Total deflection due to flexible modes at nose gear 
Slope of line drawn through two metering pin points, main gear 
Slope of line drawn through two metering pin points, nose gear 
Pneumatic force, main gear 
Pneumatic force, nose gear 
CG vertical displacement 
Unsprung mass vertical displacement, front main gear 
Unsprung mass vertical displacement, nose gear 
Rigid body pitch angle of aircraft 
Horizontal distance of aircraft 


Unsprung mass vertical displacement, rear main gear 
(Class E simulation only) 


Time derivative of T(2) 
Time derivative of T(4) 
Time derivative of T(6) 


Time derivative of T(8) 


*Class B simulation with Strut Friction only 
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SYMBOL 


TD( 10) 
TD(12) 
TDD(2) 
TDD (4) 
TOD (6) 
TDD (8) 
TOD(10) 
TD0(12) 
VELM 
VELN 
VLIFT 
XMAIN 
XMLK 
XNLK 
XNOSE 


YCEPN 


RMIT 


RM2T 


DEFINITION 

Time derivative of 1(19) 

Time derivative of T(12) (Class £ simulation only) 
Time derivative of TD(2) 

Time derivative of TD(4) 

Time derivative of TD(6) 

Time derivative of TD(d) 

Time derivative of TD(10) 

Time derivative of TD(12) (Class E simulation only) 
Total strut velocity, main gear 

Total strut velocity, nose gear 

Aerodynamic lift force 
Strut stroke, main gear 
Absolute value of XMAIN 
Absolute value of XNOSE 


Strut stroke, main gear 


Y intercept of line drawn through two metering pin points, 
nose gear 


Distance of aircraft from beginning of a 4 ft runway segment 


Runway elevation, rear main gear 


Runway elevation, front main gear (Class E Simulation only) 


Runway elevation, nose gear 


Static reaction force at main gear 


Static reaction force at rear main gear (Class E 
simulation only) 


Static reaction force at front main gear {Class E 
simulation only) 
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SYMBOL 


RN 


RSM 


RSN 


THETAI 


XRMAIN 


ZCGI 


zMl 


DEFINITION 
Static reaction force at nose gear 
RM - WM 
KN - WN 
Rigid body initial pitch angle 


Test variable for rigid body initial conditions 


(Class £ simulation only) 
Initial CG vertical displacement 
Initial tire deflection, main gear 


Initial tire deflection, nose gear 


Coefficient of polynomial fit to runway profile segment 
Coefficient of polynomial fit to runway profile segment 
Coefficient of polynomial fit to runway profile segment 
Coefficient of polynomial fit to runway profile segment 


Runway profile elevation value 


Index variable 
Number of values in metering pin - stroke table 


Metering pin diameter or net orifice area for aircraft with 
metering tubes or flutted metering pins 


Strut stroke in metering pin table 
Slape of line drawn between two metering pin points 


Y intercept of line drawn between two metering pin points 


APPENDIX E 


TYPICAL WIDE BODIED TRI-JET TRANSPORT DATA ANE SIMULATIONS 


This appendix contains data obtained late in the study which 
represents a "TYPICAL" Wide Bodied Tri-Jet Transport and may be used to 
simulate Class C and D aircraft of this type. This appendix a}so c<on- 
tains takeoff simulations over the Will Rosers and Dulles Profiles. 
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